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Abstract. We have obtained linear polarization measurements of the Large Magellanic Cloud B[e] supergiant 
S 111 using optical imaging polarimetry. The intrinsic polarization found is consistent with the presence of 
an axisymmetric circumstellar envelope. We have additionally estimated the electron density for S 111 using 
data from the literature and revisited the correlation between polarization and envelope parameters of the B[e] 
supergiant stars using more recent IR calibration color data. The data suggest that the polarization can be indeed 
explained by electron scattering. We have used Monte Carlo codes to model the continuum polarization of the 
Magellanic B[e] supergiants. The results indicate that the electron density distribution in their envelopes is closer 
to a homogeneous distribution rather than an r -2 dependence. At the same time, the data are best fitted by 
a spherical distribution with density contrast than a cylindrical distribution. The data and the model results 
support the idea of the presence of an equatorial disk and of the two-component wind model for the envelopes of 
the B[e] supergiants. Spectropolarimetry would help further our knowledge of these envelopes. 
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1. Introduction 

In the HR diagram, there is a temperature dependent 
upper limit for the stellar luminosity (Humphreys & 
Davidson 1979). Among the stars located around that re- 
gion, we find the Luminous Blue Variables (LBVs), Wolf- 
Rayet stars and the less well understood B[e] supergiants 
(B[e]SG). B[e]SG are stars with permitted and forbidden 
lines in their optical spectra dominated by strong Balmer 
emission lines, sometimes with P Cygni profiles. A char- 
acteristic of this group is a strong infrared excess that 
provides evidence for hot circumstellar dust. These stars 
show simultaneously narrow emission lines in the opti- 
cal and broad absorption lines in the UV. Such hybrid 
spectra can be explained by a two-component wind model 



a hot, fast, polar wind and a cool, slow, dense wind in the 
equatorial region. Recent reviews about B[e] stars can be 
found in Lamers et al. (1998) and Zickgraf (1999). 

At present, fifteen stars of this class are known in the 
MC and eight in the Galaxy. Data on the MC B[e]SG are 
still badly needed. For instance, B[e]SG do not show in 
general photometric variability (Am < 0.1 mag). However, 
R 4 (a B[e]SG in the SMC) exhibited a typical LBV varia- 
tion in the optical (up to 0.7 mag, Zickgraf et al. 1996a). In 
addition, the luminosities of the four new B[e]SG, discov- 
ered by Gummersbach et al. (1995) in the LMC, extend 
down to log L/Lq — 4, lower than those of the already 
known BfelSG. 



(Zickgraf et al. 1985, hereafter Z85). This wind consists of 
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Magalhaes (1992, hereafter M92) obtained polarization 
measurements of nine B[e]SG in the Magellanic Clouds 
(MC), and showed evidence for the presence of non- 
spherically symmetric envelopes, giving support for the 
two-component wind model. The available data correlated 
some with the average electron density of the envelope, N e , 
but correlated slightly better with the infrared [K — L] 
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dust excess. M92 described ways in which polarimetry of 
B[c]SG could be useful in probing their envelopes. 

HDE 269599s (= S 111 = Sk 147a, Henize 1956 and 
Sanduleak 1968, respectively) is a B[e]SG in the LMC that 
belongs to a compact cluster (~10") with about 10 stars 
(Appenzeller et al. 1984). This star could not be observed 
with photoelectric polarimetry by M92. It has however 
later been observed by us with CCD imaging polarimetry. 
In this paper, we add S 111 to the group of MC B[e]SG 
with known polarimetry (Sect. 2). After estimating the en- 
velope electron density of S 111 we revisit the correlations 
between B[c]SG polarization and N e and dust excess in the 
light of newer observational data in the literature (Sect. 3). 
We then model the intrinsic polarization of these stars by 
using Monte Carlo codes and try to infer properties of the 
B[e] envelopes (Sect. 4). Conclusions are presented in the 
last section. 

2. Observations and data reduction 

We obtained optical linear polarization images of S 111 
at the 1.5 m telescope of CTIO in December 1991. We 
used the CCD Tektronix 1024x1024 Tek#l camera and 
the measurements were made through a V filter. The po- 
larimeter is a modification of that direct CCD camera to 
allow for high precision imaging polarimetry. It consists 
basically of a rotateable half-wave plate followed by a fixed 
analyser and a filter. The arrangement is described in more 
detail by Magalhaes et al. fll996j ). 

As analyser we used a custom-made, 44 mm wide 
square double calcite block (built at Optoeletronica, Sao 
Paulo). Each component prism was cut with its optical 
axis at 45° to their faces and they were cemented with 
their optical axes crossed. This arrangement minimizes the 
astigmatism and color which are p resent when a single cal- 
cite block is used (Serkowski 1974 ). This Savart plate gives 
us two images of each object in the field, separated by 1 
mm (corresponding to about 18" at the telescope focal 
plane) and with orthogonal polarizations. By rotating the 
waveplate between CCD exposures, the ratio of the flux 
in each of the object's images changes with an amplitude 
which is proportional to the incoming beam's polarization. 
One polarization modulation cycle is covered for every 90° 
rotation of the waveplate. The simultaneous observations 
of the two images allows observing under non-photometric 
conditions at the same time that the sky polarization is 
practically cancelled (Magalhaes et al. 1996). 

CCD exposures were taken with the waveplate rotated 
through 8 positions 22? 5 apart. The exposure time at each 
position was 60 s. After bias and flatfield corrections, pho- 
tometry was performed on the images of objects in the 
field with IRAF, and then a special purpose FORTRAN 
routine processed these data files and calculated the nor- 
malized linear polarization from a least squares solution. 
This yields the Stokes parameters Q and U as well as the 
theoretical (i.e., photon noise) and measurement errors. 
The latter are obtained from the residuals of the observa- 
tions at each waveplate position angle (^j) with regards to 
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Fig. 1. The HDE 269995 cluster in the LMC, from a sec- 
tion of an image taken in CTIO. S 111 is star number 
1. The others were used for determining the foreground 
interstellar polarization towards S 111. 



the expected cos Aipi curve and are quoted in Table 1; they 
are consistent with the photon noise errors (Magalhaes et 
al. ~ 



1984). The instrumental Q and U values were con- 



verted to the equatorial system from standard star data 
obtained in the same night. The instrumental polarization 
was measured to be less than 0.03%. 

To obtain the intrinsic polarization we need to know 
the foreground interstellar polarization. In order to do 
this, we employed field stars. The sample that we used as 
field stars corresponds to three stars in the same cluster 
as S 111 (Fig. 1). The advantage of using them is because 
we know that those stars are at the same distance from 
us as S 111 as well as angularly close to it. The intrinsic 
polarization is calculated as follows: 



Qi — Qo — Qis, 
Ui = U - U m 



(1) 



where the (Qis, £7is) values are weighted mean of the three 
stars. Subscripts i, o and IS correspond to intrinsic, ob- 
served and interstellar polarization, respectively. 

In Table 1 we show the observed polarization measure- 
ments (P — \J Q 2 + U 2 ) for four stars of the cluster. The 
stars labeled with 2, 3, 4 in Fig. 1 were used as field stars 
to calculate the foreground interstellar polarization. Sill 
is star 1. In Table 2 we present the optical polarization 
data for S 111. The last 3 columns of Table 2 show that 
S 111 has intrinsic optical polarization. 

3. Polarization and the B[e]SG envelopes 

The envelopes of B[e]SG provide an ample opportunity 
for scattering of radiation from the central star due to 
the available free electrons, as revealed from their emis- 
sion lines and near IR (J band) excess, and dust, evi- 
denced by their IR excess. Following M92, we reinvesti- 
gate the correlations between the optical polarization of 
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Table 1. Observed Optical Polarization Data for the 
Cluster 



N 


P(%) 


OP (%) 


e n 


1 (=S 111) 


0.710 


0.033 


50.5 


2 


0.672 


0.085 


26.1 


3 


0.532 


0.075 


35.8 


4 


0.544 


0.029 


29.7 



B[c]SG and inferred electron densities and infrared ex- 
cesses. This reanalysis is warranted by the additional data 
point Sill provides, the recent finding concerning the bi- 
narity of AV 16/R 4 and newer photometric calibration 
tables for optical-near IR photometry. 

3.1. Polarization and Electron density 

It is believed that the J excess from B[c]SG is due to free- 
free and bound-free radiation. Therefore, such excess is a 
measure of the electron density in the envelopes of B[e]SG 
(Zickgraf et al. |1986| , hereafter Z86). M92 correlated the 
polarization of nine MC B[e]SG with the electron density 
(N e ) obtained from J excesses (AJ) by Z86. A measure of 
the optical magnitude needed to fit the stellar continuum 
energy distribution, and hence for obtaining AJ, is not 
available for S 111. Thus we estimated N e for S 111 not 
from A J but from the emission measure (EM) of Ha. 

For estimating N e from Ha we followed the method 
outlined by Z86. As a check on our procedure, we initially 
applied it to the other B[e]SGs, for which optical mag- 
nitudes are available. The first step required is obtaining 
the Ha continuum flux. This was done using an interpola- 
tion of the fluxes in the V and / filters using the observed 
magnitudes of Z85, Z86 and Zickgraf et al. 1989 (hereafter 
Z89) and a power law, F\ = fcA~ n . The exponent n was 
found to vary between 1.5 and 2.9 (with median = 2.5) for 
the different stars. The magnitudes were corrected for in- 
terstellar extinction, using the extinction curve of Code et 
al. (1976). The Galactic color excess E(B — V) was t aken 
as 0.05 (Z85, Bessell 1991 and Rodrigues et al. |1997| ) and 
the Magellanic color excess E(B — V) from Z86 and Z89. 
The fluxes were then computed using the absolute fluxes 
of Bessell (1979). Our results compared well to those of 
Z89 within 30%. 

For S 111, in order to calculate the Ha continuum flux 
we estimated the magnitude of the object in the V fil- 
ter, using the bolometric correction BCv and the stellar 



luminosity by McGregor et al. ( |1988| , hereafter M88), as 
well as the distance modulus for the LMC (Walker 1998). 
The estimated intrinsic (i.e., extinction corrected) appar- 
ent magnitude for S 111 thus found was 9.20 mag. This 
estimated V magnitude is not quite consistent with the 
near infrared magnitudes of M88 (their Table 5). We did 
perform another estimate of the magnitude of S 111 us- 
ing the photometry of Mendoza (1970), who obtained the 
magnitude in the V filter for the cluster to which S 111 
belongs. Our CCD images provide the relative fluxes of 
the brightest stars in the cluster. It was then possible to 
estimate the magnitude of S 111 from the cluster V mag- 
nitude of Mendoza (1970) and our relative flux data. For 
this purpose, we assumed that the flux of the cluster was 
due to the contribution of its four brightest stars. From 
this assumption we obtained a new intrinsic magnitude for 
S 111, 10.26 (with an error of ~0.02, the error in V for the 
Sill cluster as estimated by Mendoza), a value more con- 
sistent with the near infrared magnitudes of S 111. Since 
there is no available / filter measurement for Sill either, 
the Ha continuum flux of S 111 was obtained from an ex- 
trapolation of this V flux assuming a value of n=2.5 for 
the index of the continuum flux distribution. A more def- 
inite estimate should probably be best done from direct 
filter photometry for S 111, specially in the I filter. 

The next step was obtaining the Ha flux (Ha contin- 
uum flux x equivalent width) and then the Ha luminosity. 
For S 111, the equivalent width was measured from the Ha 
profile given by Stahl & Wolf (1986) as 356 A. The EM 
was calculated using the case B recombination theory for- 
mula (Osterbrock 1989 ). From the definition of EM, N e 
can be calculated as -\/EM/V, where V is the volume of 
the envelope. The shape of the envelope was taken as sug- 
gested by Z85. It consists of a disk with inner radius of 
1 i?*, outer radius of 300 R± and thickness of 1 i?*. For the 
stellar radius of S 111, we obtained the value of 148 Rq us- 
ing the effective temperature and stellar luminosity from 
M88. The resultant N e for S 111 was 0.43 x 10 9 cm" 3 . 

As a check on the procedure, we reevaluated N e for 
other B[e]SG. The Ha equivalent width of other B[e]SG 
stars were either taken from the literature (R 126, S 134 
and S 12; Z89) or estimated by us (S 18, S 22 and R 66). 
Among the latter, two compared very well with the litera- 
ture (S 18, Z89; S 22, Schulte-Ladbeck and Clayton 1993). 
Both the fluxes and EM values obtained by us agreed with 
those of Z86 and Z89 estimated from Ha. On the other 
hand, our values of N e differ some, though always of the 
same order of magnitude, from those of Z89's obtained 
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from J excess. In general, the electron density from AJ 
is on average two times larger than that of Ha. We hence 
used 0.8 x 10 9 cm~ 3 for the electron density in the enve- 
lope of S 111. 

As discussed in detail by M92, high levels of intrin- 
sic polarization are related to edge-on stars. We hence 
expect that polarization of these stars might be propor- 
tional to the electron density and the other stars (pole-on 
and intermcdiatc-on) be distributed in a triangular pat- 
tern due to the dependence on sin 2 i, where i is the incli- 
nation angle through which the discs are seen. In Fig. 2 
we plot the intrinsic polarization as a function of electron 
density. We use our N e value above for S 111 and Z89's 
for the remaining stars. The polarization values are taken 
from Sect. 2 of this paper for S 111 and from M92 for the 
remaining objects. 

Unlike M92, we have omitted from the plot in Fig. 
2 the SMC star AV 16/R 4. Zickgraf et al. (1996a) have 
shown that this object is in fact a binary star, consisting of 
a B[e] star and an early A- type companion. The observed 
polarization of an object is basically the ratio between the 
scattered, polarized light from the circunstellar material 
and the total light from the system. Hence, such source 
of additional, unpolarized light in the system (i.e., the A 
component) would act to dilute an otherwise larger in- 
trinsic polarization. In fact, of the four equator-on cases 
(as evidenced by their spectroscopic properties from the 
work of Zickgraf and collaborators) analysed by M92, R 4 
showed the smallest intrinsic polarization. Interestingly, 
in addition to that, in the P vs. N e correlation R 4 of 
M92 (his Fig. 4) R 4 clearly stood out from the other 
three equator-on cases by being well off the line defined 
by those objects. 

We did a linear % 2 fit for the remaining spectroscop- 
ically edge-on stars (R 50, R 82 and S 22, filled squares 
in Fig. 2) and present that fit in Fig. 2. It can be seen 
that the remaining stars, i.e., the pole-on and intermedi- 
ate stars, are distributed below this fit line, as expected. 
As in M92, we chose to use the Spearman rank-order cor- 
relation coefficient (r s ), to which a statistical significance 
can be attached. The calculated value of the r s for the P 
versus N e correlation for all objects is 0.55, with a 12% 
probability, where a small probability means high signif- 
icance. This indicates a moderate correlation between P 
and N e . Nevertheless this correlation is much higher than 
that of M92 (0.32 with 41%), which included R 4. 

3.2. Color excess A(K — L) 

We can use the near-infrared color excess at wavelenghts 
longer than J as a measure of dust emission (Z86, Z89). 
We define the following color excess A(K — L): 



A(K -L) = (K- L) corr - (K - L)i 



(2) 



a R66/S73/Sk-6946 
b R50/S65/Sk193 
c S 18 

d R126/S127/Sk-69 216 
e S 12 /Sk -67 23 
f S 134 /Sk -69 259 
g R82/S89/Sk-69 77 
h S22/Sk-67 64 



i S 111 



where (K — L) COII is the observed color index (K — L) 
corrected for interstellar extinction and (K — L)j nt is the 
intrinsic color index due to the stellar contribution. 



N e (x 10 cm ) 

Fig. 2. Intrinsic polarization versus electron density. The 
dotted line corresponds to a \ 2 linear fit for edge-on stars 
(R 50, R 82 and S 22), indicated by the filled squares. 



(K — i)i„t was taken from the newly available 
UBVRIJHKL color tables of Bessell et al. (1998). These 
tables take into account the star's gravity, finer intervals 
of effective temperature and colors for T > 21000K com- 
pared to those of Johnson (1966). (K — L) COII was obtained 
using the observed magnitudes of Z85, Z86 and Z89, cor- 
rected by the interstellar color excesses taken from Z85 
for the Galaxy and from Z86 and Z89 for the MC and us- 
ing the extinction coefficients of Leitherer & Wolf (1984). 
We note that there might be probably uncertainties in the 
color excesses derived in these papers. The reason is that 
the dust envelope may also contribute to the observed 
reddening of the stellar spectrum. The amount of such 
reddening will depend on detailed modeling. The derived 
interstellar color excesses may then be considered as upper 
limits. One such example may be S 18, which has a derived 
E(B - V) = 0.40 (Z89). Such a relatively large color ex- 
cess, albeit possible, is rarely observed towards the SMC. 
Nevertheless, as expected, these color excesses have only 
a small effect in (K — L) COII at these long, IR wavelengths. 

In Fig. 3 we plot the intrinsic polarization as a function 
of the color excess A{K — L). As before, we omit R 4 from 
the correlation. A linear \ 2 n t to the spectroscopically 
edge-on stars (R 50, R 82 and S 22) is also shown in the 
figure. It can be seen that not all stars are distributed 
below the fit line. The Spearman correlation coefficient 
for the P versus A(K — L) relation for all objects is 0.20, 
with a 60% probability. This value of r s shows a very weak 
correlation between P and A(K — L). 

In summary, the existing B[c]SG polarimetric data 
show that intrinsic polarization correlates better with elec- 
tron density than dust scattering. This result is consis- 
tent with the spectropolarimctric observations of S 22 by 
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a R66/S73/Sk-69 46 
b R50/S65/Sk193 
c S 18 

d R126/S127/Sk-69 216 
e S 12 /Sk -67 23 
f S 134 /Sk -69 259 
g R82/S89/Sk-69 77 
h S22/Sk-67 64 



J S 111 



Table 3. Electron Densities and Stellar Radii for MC 
B[e]SG 



1.6 1.1 

Delta (K-L) (mag) 



2.0 



Fig. 3. Intrinsic polarization versus color excess. The dot- 
ted line as in Fig. 2. 

Schulte-Ladbeck & Clayton (1993). They found that the 
polarizing mechanism in that object can be explained by 
electron scattering. Additional multifhter and spectropo- 
larimetric measurements of the B[e]SG are clearly of the 
greatest interest. As discussed in M92, these measure- 
ments will allow an improved separation between intrinsic 
and foreground polarization as well as the detection of 
changes in polarization across spectral features produced 
in the B[c]SG envelopes. Such changes will help reveal the 
physical structure of their envelopes. 

4. Monte Carlo models 

In order to quantitatively infer consequences of the ob- 
served polarization to the B[e]SG envelopes, we used a 
Monte Carlo code developed by our group (Carciofi & 
Magalhaes 2001). In this preliminary study, two geome- 
tries were chosen for these envelopes (although others 
may be used): cylindrical geometry and spherical geom- 
etry with density contrast. The density laws used were 
one of constant density and one proportional to r~ 2 . The 
codes perform the radiative transfer of the Stokes param- 
eters (I, Q, U, V) for photons in scattering envelopes. 

From the preceeding section, we found that polariza- 
tion correlates best with electron density. We have hence 
only considered models with electron scattering at this 
point. We used data on nine B[e]SG for the models; that 
is, those from M92, with the exception of the binary 
AV 16/R 4, plus S 111. In Table 3 we summarize their 
stellar radius (Zickgraf 1990) and electron density values 
(Z89). For the modelling we considered five values for the 
average electron density in the envelope which encompass 
the observed ones (Table 3): [0.5, 1.5, 3, 4, 6] x 10 9 cm" 3 . 
The mean electron density for each case is calculated from 
< jV e 2 > (= EM/V), as described in Sect. 3.1. For the mod- 
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S 111 
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R 126 


72 
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S 134 


45 
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S 12 


30 


3.3 


R 66 


125 


1.2 


R 82* 


50 


2.8 


S 22* 


49 


6.0 


S 18 


35 


3.4 


R 50* 


81 


1.5 



a stellar radii (Zickgraf 1990) 

b electron density values (Z89) 

c values obtained in this work 

* edge-on stars 



els, we have used the same approach, i.e., we calculated 
< iV 2 > for a given geometry and density law and then 
took the square root of this value to represent the mean 
electron density of that specific model. In the figures of 

1/2 

the models we use < iV 2 > along the x-axis. 
4.1. Cylindrical geometry 

The envelope with cylindrical geometry was modeled 
adopting a stellar radius of 70 Rq, typical of a B[e]SG. 
The dimensions of the cylindrical disk are the same used 
earlier (Sect. 3.1), namely, inner radius R4 = 1 i?*, outer 
radius R = 300 i?* and thickness of 1 i?*. 
The electron density law used was: 



n e (r) 



K 



(3) 



where K is a constant and n is an exponent for constant 
(n=0) or variable (n=2) densities across the envelope. The 
constant K can be derived using the relation between the 
optical depth and the density law, 



Ro 

t = J <j^n e (r)dr . 

Ri 



(4) 



with <tt being the Thomson scattering cross section. 
4.1.1. Homogeneous density 

In Fig. 4 we show the intrinsic polarization as a function 
of mean electron density (and equatorial optical depth) 
for a cylindrical disk with homogeneous density for sev- 
eral inclination angles. From the figure, we can note that 
these models fail to fit the observational data as most data 
points, including the edge-on objects, lie well above the 
curve i=90°, which is the maximum polarization predicted 
by the model. We conclude that homogeneous cylindrical 
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a R66/S73/Sk-6946 
b R50/S65/Sk193 
c S 18 

d R126/S127/Sk-69 216 
e S 12/Sk-67 23 
f S134/Sk -69 259 
g R82/S89/Sk-69 77 
h S22/Sk-67 64 



S 111 ! 



2 (x10 9 cm 3 



Fig. 4. Intrinsic polarization of MC B[e]SG as a function 

1 /2 

of envelope mean electron density N e (= < N 2 > ) and 
equatorial optical depth r. Lines represent models with 
cylindrical geometry and homogeneous electron density 
for several inclination angles. The points with error bars 
correspond to polarimetric data from this paper (S 111) 
and M92 (other objects). 

envelopes are not a good representation for the MC B[e]SG 
envelopes. 

4.1.2. Density proportional to r~ 2 

The model results when the density falls with r~ 2 are 
presented in Fig. 5. In the figure, most of the observa- 
tional data points lie between the inclination angle curves 
i=20°and i=45°. In particular, the edge-on objects fall far 
from the i=90°curve. We conclude that, like the constant 
density models, the cylindrical r~ 2 models do not fit well 
the observations either. 

In order to investigate the dependence of the model 
results on the assumed stellar radius, we ran two addi- 
tional cases, i?*=30 Rq and i?+=140 Rq, which reason- 
ably cover the range of MC B[e]SG radii (Table 3), for 
both constant and r~ 2 density models. For the cases of 
constant density with i?*=30 Rq, the model curves fall 
below all observed data points. The i?*=140 Rq model 
curves cover the P — N e data space better but still fail to 
adequately fit the edge-on cases. In addition, we know the 
edge-on objects do have in fact i?+ between about 50 and 
80 Rq (Table 3). For the r~ 2 density models, neither of 
these radius values fit the data points adequately. 

In summary, we conclude that the models using cylin- 
drical geometry do not fit with the data for B[e]SG, at 
least not for the dimensions that have been proposed 
(Z85). 
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Fig. 5. Intrinsic polarization of MC B[e]SG as a function 
of mean electron density square root < iV 2 > and equa- 
torial optical depth r. Data are as in Fig. 4. Lines repre- 
sent models with cylindrical geometry and electron density 
varying with r~ 2 for several inclination angles. 



4.2. Spherical geometry with density contrast 

To model the spherical geometry with density contrast, we 
chose to use the parameterization of Waters et al. (1987), 
which has been also used by Wood et al. (1996). The 
equatorial-to-polar density ratio may in fact render this 
model more realistic for the B[e]SG. 

The Waters et al. (1987) density law is 



( ) i:i + Lsii."'fl). 



(5) 



where no is the density in the stellar surface, 9 is the polar 
angle. This parameterization consists of a spherically sym- 
metric region combined with an equatorial disk in which 
(1 + A) determines the density contrast, between the pole 
and equator, and m is a measure of the opening angle of 
the disk. The aperture angle at which the density is half 
its maximum value (measured from the equator) is: 



0o 



cos 



-1/Iu/m 



(6) 



We have chosen to model two cases, that of an enve- 
lope with a constant density (n=0), and that with density 
proportional to r~ 2 (n=2). In these cases, we used the 
same amount of mass as that of the scattering cylindrical 
disk. In doing that, the spherical envelope outer radius 
derived was 2850 Rq (~41 R±) . Again, the stellar radius 
used in the models was 70 Rq. 
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a R66 
b R50 
c S 18 
d R 126 
e S 12 
f S 134 
g R82 
h S22 



hi i = 90° 



,i i = 60° 



2 4 6 

<N2> 1/2 (x 10 9 cm 3 ) 

Fig. 6. Intrinsic polarization of MC B[e]SG as a function 

1 /2 

of electron density N e (— < N 2 > ) and equatorial op- 
tical depth t. Data point are as in Fig. 4. Lines represent 
models with spherical geometry with density contrast and 
uniform electron density. 




Fig. 7. Intrinsic polarization of MC B[e]SG as a function 
of mean electron density square root ^J<N%> and the 
equatorial optical depth r. Data are the same of Fig. 4. 
Here we considered a spherical geometry with density con- 
trast, and a r~ 2 electron density. 



4.2.1. Constant density 

In the case of homogeneous density, we see that in eq. 5, 
n e (r) only depends on 9, m and A. We set the opening an- 
gle at 0=5°, for which ra=182. Zickgraf (1989) estimated 
disk opening angles of 20° to 40°. We have run models 
with this range of values but the resulting polarization 
leves were too high (e.g., 7% for N e = 6 x 10 9 cm -3 ). 
The value of density contrast used was 10 3 (Z86), which 
implies in A=999 in eq. 5. 

In Fig. 6, we present the results for envelopes with the 
geometry parameterized with eq. 5 and constant electron 
density. These models provide a good fit to the observa- 
tions, including the edge-on objects. We note that polar- 
ization level for i=90° is larger than that obtained with 
cylindrical geometry (Fig. 4). This is due to the star, in 
the cylindrical geometry case, not being completely oc- 
culted by the envelope. The star then contributes with 
direct, non-polarized flux, lowering the resulting polariza- 
tion level. 

4.2.2. Density proportional to r~ 2 

With the density now proportional to r -2 , we use the same 
parameter set as for constant density. The model results 
are presented in Fig. 7. In the figure we can see that the 
model polarization obtained fall well below the expected 
level. 



5. Conclusions 

We present optical polarization data for the LMC B[e]SG 
S 111. The estimated intrinsic polarization is (0.467 ± 
0.042)% at 75?3. This polarization is consistent with the 
presence of a non-spherical circumstellar envelope around 
the object. 

After excluding AV 16/R 4, a B[e]SG now known to be 
a binary, and with the inclusion of S 111, we found using 
more recent IR calibration data that intrinsic polarization 
correlates better with the electron density rather than the 
color excess A(K — L) in the envelopes of B[e]SG. This 
suggests that electron scattering is the main mechanism 
that causes polarization in the envelopes of these stars. 

The spectroscopic data have been suggestive (Z85) of a 
two-component wind model for B[e]SG, i.e., a slow equato- 
rial wind and a fast polar wind. We have hence considered 
two idealized models for such wind, one with a cylindrical 
geometry and another where the material is spherically 
distributed around the star but with a contrast enhance- 
ment along the equator. The results of our Monte Carlo 
modeling show that the spherical geometry with such den- 
sity contrast fits the polarization observations very reason- 
ably, being consistent with the two-component model for 
B[e]SG. In addition, the constant density case fits the data 
better than an r~ 2 electron density law. This result agrees 
with the presence of a slow equatorial wind observed and 
modeled by Zickgraf et al. (1996b). 

Questions still to be answered include details of the 
geometry of the equatorial disk, such as its opening an- 
gle, and the precise location of the scattering material 
within the disk. These will have to await spectropolari- 
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metric observations of these objects. Changes (or the ab- 
sence of them) in polarization accross spectral features, 
e.g., emission lines, would help us point out more pre- 
cisely where the polarization arises in the envelope (de- 
creases across emission line, for instance, would indicate 
that the polarization is produced closer to the star, etc.). 
Other spectropolarimetric features, such as changes across 
the Balmer discontinuities, might be helpful in assessing 
the envelope's opening angle (e.g., Wood et al. 1997, for 
Be stars). 
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